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1 Introduction 


In flow ducts, transient acoustic waves are rapidly dissipated by viscous forces that are 
present in all real flows. This phenomenon leaves the periodic sound as the dominate source 
of noise for most duct acoustics problems. Understanding the behavior of these periodic 
sound waves has long been an important goal in duct acoustics research. The importance 
of this goal has further increased as aircraft acoustic analyses have been required to include 
the effects of increasingly severe operating conditions. To date, most analytical research has 
used simplified frequency-domain models to obtain knowledge in regard to flow effects on 
periodic sound waves. However, the research community has generally accepted the fact that 
a fully adequate treatment of sound effects in flow ducts will be achieved only through a 
direct numerical solution to Euler equations. In this approach, periodic acoustic waves are 
obtained by integrating the Euler equations to a periodic state. Unfortunately, the absence 
of dissipation from these equations means that acoustic transients that are introduced when 
these fluid dynamics equations are integrated from a given initial state do not decay naturally 
as they do in reality. As a result, transient waves must be eliminated by other means so that 
the proper periodic state can be achieved. 

One means for eliminating initial transients from a set of inviscid equations such as 
the Euler equations is to include artificial viscosity in the numerical scheme. However, 
the inclusion of the artificial viscosity often damps the acoustic disturbances so that the 
appropiate periodic steady state is not captured. The method often used to extract the 
periodic steady-state field from Euler’s equations is to specify time-domain conditions on the 
boundary of the computational domain that guarantee that incident energy is not reflected. 
This subject was treated, for example, in references [l]-[6]. The computational boundary 
conditions that axe presented in these works are local; and these boundary conditions are not 
generalizable to higher dimensions, arbitrary sound sources, and acoustically treated ducts. 

In two recent papers (references [7, 8]), nonreflecting boundary conditions were shown to 
be obtainable in the frequency-domain by using nonlocal computational boundary conditions. 
The purpose of this paper is to extend these frequency-domain computational boundary 
conditions to the time domain. When nonreflecting frequency domain boundary conditions 
are applied in the time domain, they are not necessarily transparent to transient waves 
because, in time-domain computations, transient acoustic waves are introduced by the initial 
conditions. Thus, one issue that must be dealt with immediately in this context is the 
elimination of transient acoustic waves, which were absent in the frequency-domain analysis 
of the two earlier works (reference [7, 8]). Here, a set of time-domain computational boundary 
conditions that allow periodic waves to pass without reflecting is derived. These boundary 
conditions are referred to as periodic nonreflecting time-domain boundary conditions. 

In the following section the physical duct acoustics problem, which is infinite in ex- 
tent, is described and the initial boundary value problem is formulated. In section 3, the 
time-domain computational boundary-condition operators are presented. These boundary- 
condition operators allow the domain to be truncated to a finite length to achieve computa- 
tional efficiency. These conditions are shown to be convolution integrals of the corresponding 
frequency- domain operators that are presented in reference [7]. In section 4, the boundary 
conditions are coupled to the Euler equations and advanced to a periodic state by using 
MacCormack’s algorithm. Results for both a pure and a multitone source for several initial 
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conditions axe given in section 5. Conclusions are given in section 6. 


2 Problem Formulation 

2.1 Rectangular Duct and Coordinate System 

Consider a rectangular duct with the inflow and outflow boundary at x = 0 and x = L, 
respectively, as shown in figure 1. Although the duct is assumed to be infinitely long in the 
axial direction, the computational domain is limited to a finite length L, and the sound field 
is considered to be independent of the y direction (i.e., planar). The sound source is located 
at the inflow boundary and is assumed to be periodic in time with a fundamental period of, 
T. The sound source generates acoustic disturbances that propagate from left to right in the 
duct. Only compressible, inviscid disturbances without heat conduction are considered. 

For purposes of this paper, the fluid in the duct is assumed to decompose into a back- 
ground and an acoustic component. The background pressure Pq, density po, and tempera- 
ture To are constant. The background flow has zero velocity. The cross-sectional area of the 
duct is constant, although extension of the analysis to higher dimensions and variable-area 
ducts is straightforward. At time t = 0, the disturbance in the duct is assumed to be known. 
As is common in acoustics, the source pressure g(t ) is assumed to be known for all t > 0 and 
drives the flow. Finally, acoustic waves that exit the computational domain at x = L do so 
without reflecting. The goal here is to derive time-domain nonreflecting boundary conditions 
that allow the duct to be terminated at x = L. These boundary conditions must also be 
capable of driving the system to a time-periodic state. 


2.2 Governing Equations 

The fundamental equations that govern inviscid compressible disturbances in the duct are 
the Euler equations (reference [9]): 


m } dm 
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where C v is the specific heat at constant volume. Equation (1) is a vector equation that 
contains three components. Physically, the first, second, and third components represent the 
conservation of the mass, the x momentum, and the energy, respectively, within the duct. 
The unknown variables are the fluid density p, axial velocity u, temperature T, and pressure 
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p. Thus, the vector equation (1) is a system of three equations in four unknown variables 
p,p,T, and u. To close the system, an ideal gas is assumed: 

P = pRT (5) 

Here, R is the universal gas constant. Equation (5) is used to eliminate temperature and to 
reduce the system to three equations in three unknown variables. 

Equations (l)-(5) are the most general set of nonlinear fluid dynamic equations; these 
equations govern inviscid compressible disturbances in the duct. No known periodic solutions 
exist for these equations in this most general form. Therefore, the objective in direct numer- 
ical simulation is to numerically evolve these equations to a periodic state which requires the 
specification of initial and boundary conditions. 

2.3 Initial and Boundary Conditions 

As is common in initial boundary value problems in fluid dynamics, the flow will be initiated 
from a state of rest: 

p(x, 0) = P 0 (6) 


p(x,0) = p o (7) 

and 

u(r,0) = 0 (8) 

where Pq and p 0 are constants. The governing equations require inflow and outflow boundary 
conditions. The Euler equations in one space dimension and time require two boundary 
conditions at the subsonic inflow and one condition at the subsonic outflow (reference [3]). 
At the inflow, the disturbance pressure is prescribed so that 

P(0, t) = Po + g(t) (9) 


The inflow acoustic disturbance is assumed to be linear and homentropic: 

dp(0,t) 2 dp(0,t) 

~Si = ° 

Here, cq is the ambient speed of sound which is defined as 


( 10 ) 



( 11 ) 


where 7 is the ratio of specific heats and is approximately equal to 1.4 for air at standard 
conditions. Finally, the initial boundary value problem is fully specified when a radiation 
condition is specified at the outflow boundary x = L. Development of the radiation condition 
is the subject of the following section. 
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3 Radiation Condition for Periodic Solutions 


The radiation condition relates the acoustic pressure p and the normal component of acoustic 
velocity u at the outflow boundary x = L. The Fourier series for the periodic pressure and 
axial velocity can be given in complex form as 


p(L,t) = P 0 + p{L,t) 

(12) 

u(L,t ) = Uo + u(L,t) 

(13) 

OO 

P{L,t) = J2 p(£;w„)e ,u '" t 

n=— co 

(14) 

OO 

u(L,t)= ^2 u(L]u> n )e ,Wn< 

(15) 


n~ — OO 


the frequency u n is defined in terms of harmonic number n 

and the fundamental period of 

the sound source T: 

2nx 

“n = j, 

(16) 

where n is an integer. The periodicity condition is 


p(L,t + T ) = p(L,t) 

(17) 

and 


u(L, t + T) = u(L , t) 

(18) 

Band-limited periodic functions contain a finite number of terms n < N such that 

P(L,t) = J2 p(L-,u> n )e' Wnt 

(19) 


n=-A r 


«(£-,<)= ^2 u(L;u> n )e tUnt ( 20 ) 

n=—N 

The pressure and axial velocity used here must each be physically real-valued, so that the 
series coefficients satisfy the conjugate identity for negative frequencies: 

p(L-,-u; n )=p*(L; u n ) (21) 

and 

u(L\ -w n ) = u*(Z,;u>„) (22) 
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where the superscript * denotes the complex conjugate. These coefficients are integrals (over 
a single period) of the product of time functions and the complex exponential e~ lwt : 

P{L ; oj n ) = ^ e~'“ nT p{L, t ) dr (23) 

and 

u(L] w n ) = ^ e~ tUnT u(L, r) dr (24) 

In this definition of the coefficient, the choice of t is unimportant; however, the definition has 
the physical implication that the coefficient depends on the time period that immediately 
preceeds the present t. The band-limited pressure and normal velocity functions considered 
here are not restricted to zero mean values, so that the zero-frequency coefficient is not 
necessarily zero: 


p{L\u 0 ) ^ 0 

(25) 

u(L]u> 0 ) ^ 0 

(26) 


This requirement is necessary to accommodate nonlinear effects. Note that only ( N + 1) 
complex constants are needed to represent each fluctuating function p(L,t ) and u(L,t). 
Often, the integer N is taken to be some integer power M of 2 to take advantage of the 
FFT: 

N = 2 m (27) 

In practical problems, M may range from 5 to 10, so that N would range from 32 to 1,024. 
(These numbers are only typical values.) 

The boundary conditions utilized here use generalized forms of impedance and admit- 
tance. The elementary scalar forms of these relations in the frequency-domain are 

p(L;uj n ) = Z(u n )u(L]u> n ) (28) 

u(Z;u>„) = (3(u n )p(L;u > n ) (29) 

where Z(u > n ) is the impedance and /3(o>„) is the admittance. Values of Z{u n ) and /3(a; n ) 
that correspond to outgoing acoustic waves are obtained from standard frequency-domain 
analyses (reference [7]) The time- domain forms of equations (28) and (29) are given by 
convolution integrals: 


p(L,t) = 

f Z(t — t)u(L, T)dr = < Z *u> 
Jt-T 

(30) 

II 

*3 

A 

* 

V 

II 

T 

l 

hi 

(31) 

The time-domain, or temporal, impedance and admittance are functions of the time delay 
t — t, as in equation (19): 

Z(t- 

')4 e 

1 n——N 

(32) 
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( 33 ) 


W - r) = I £ 

1 n=-N 

Either equation (30) or (31) is a time- domain radiation condition to be applied at x = L. 
Results presented in this paper were obtained with the impedance form of the radiation condi- 
tion (equation (30)). In time- domain computations, the impedance or admittance function, 
defined by either equation (32) or (33) is tabulated and held for use by the computational 
algorithm. If the computation proceeds in uniform time increments, direct table lookups 
supply the requested values at minimal cost. We assume here that the preliminary com- 
putation produces a table with acceptably fine resolution, such that direct lookup supplies 
the needed value of the impedance or admittance. Uneven time increments involve addi- 
tional computation at each time step, evaluation of the series at each point, or additional 
preliminary computations to produce an enlarged preliminary tabulation that is suitable for 
interpolation by some simple formula. 

A similar argument holds for the acoustic pressure or axial velocity required to perform 
the time integration in equations (30) and (31). One of these variables must be computed by 
using the computational algorithm and then held for use for one period; the other is computed 
from radiation condition equation (30) or (31). Because the integrand of either radiation 
condition is known only at discrete time increments, numerical integration is required to 
implement the boundary condition. The advantage of the boundary-conditionformulation 
that is presented here is the generality that it affords. The convolution equations provide a 
direct means for extending the time-domain boundary condition in one space dimension to 
any number of spatial dimensions and to ducts with acoustically treated walls and arbitrary 
sources. This extension requires the use of the frequency-domain nonlocal nodal impedance 
matrix operators. These matrix impedance operators have been shown to be an accurate 
and effective computational boundary condition in duct acoustics (references [7, 8]). 

4 The Numerical Method 

Now that we have specified initial and boundary conditions for the Euler system (equa- 
tion (1)), we focus on a numerical solution technique. One specific algorithm that has 
gained wide acceptance for solving time-dependent problems in fluid dynamics is MacCor- 
mack’s predictor-corrector scheme (reference [9]). This algorithm is used here to integrate 
the Euler system, coupled with the initial, inflow, and radiation conditions to a periodic 
state. We describe the method here for the interested reader. 

To apply MacCormack’s method to the Euler system, the continuum region of interest 
(x,<), is replaced with a discrete finite-difference mesh (xj, t n ) as shown in figure 2. Thus, in 
the finite-difference mesh, we denote the pressure at grid point (ij, t n ) by p hn . The finite- 
difference mesh is assumed to contain NX evenly spaced points in the axial direction, as 
shown. The grid spacing in the x and t direction is denoted by Ax and At, respectively. 
The grid spacing in the t direction (i.e., At) is not necessarily evenly spaced. MacCormack’s 
solution of equation (1) within the grid in figure 2 takes the form of a predictor-corrector 
technique by using forward differences on the predictor step and backward differences on the 
corrector step. Although the differences are first-order accurate in each step, the overall solu- 
tion is second-order accurate in both space and time. Because MacCormack’s time-marching 
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algorithm is explicit, roundoff error can grow in an unbounded manner and destroy the solu- 
tion if the Cauchy-Friedrichs-Lewy (CFL) criterion is not satisfied. The physical significance 
of the CFL criterion is that the time step At must be, at most, equal to the time required 
for a sound wave to propagate between two adjacent grid points. The criterion gives the 
following time step restriction (reference [9]): 

At = <r t (At) max (34) 

where (A<) max is the maximum value of within the computational domain, and a t is a 
positive constant less than unity. 

The following is a summary of the steps involved to advance the solution forward in time. 

1. Choose a coordinate grid Ax that is sufficiently fine to resolve all spatial gradients in 
the duct. 


2. Beginning with the initial condition, apply equation (34) at each grid point to deter- 
mine a minimum time step At for the integration. 

3. Beginning with the known solution at t, integrate the Euler system (equation (1)) by 
using MacCormack’s scheme at all interior points to advance the interior solution to 
the next time step t + At. 


4. Special care is required to advance the solution in time along the inflow and outflow 
boundary. Along these boundaries, some conditions are prescribed; others must be 
obtained from the compatibility conditions of the Euler system. Along the inflow 
boundary x = 0, the two inflow conditions and the axial momentum equation are used 
to advanced the solution in time. To compute the inflow pressure, we use the discrete 
form of the known source condition (equation (9): 


Pi,n+i — Po + g n + i (35) 

The homoentropic condition (equation (10)) now gives the inflow density from the 
known interior solution and the inflow pressure obtained from (equation (35)) 


Pl,n+1 ~ />2,n+l 


+ 


Pl.n+1 — P2,n+1 
r 2 

°0 


(36) 


To advance the inflow axial velocity to the next time step, we integrate the axial 
momentum equation (i.e., the second component of equation (1)) in time along the 
inflow boundary with MacCormack’s scheme. This integration is performed by using 
forward differences for all axial gradients in the equation. Along the outflow boundary, 
the density p and axial velocity u are advanced in time by integrating the first two 
components of equation (1) forward in time with MacCormack’s scheme. Backward 
differences are used to approximate all axial gradients in the integration. The radiation 
condition (equation (30)) provides the necessary condition for obtaining the outflow 
pressure p: 

PNX, n+1 = Po+ < Z * UffX > (37) 

The time integration in (equation 37)) is performed numerically. Axial velocity values 
are held by the algorithm for one time period, and cubic spline interpolation is applied 
to perform this numerical integration. 
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5. Continue to advance the solution in time until a periodic state is achieved. 


5 Results and Discussion 

A computer code that implements the boundary condition derived in this paper has been 
developed and programmed to run on a DEC-Alpha work station and a Cray computer. 
Results presented in this section were obtained with the underlying objective of comparing 
the periodic waveform obtained from the boundary condition with the exact solution for a 
range of source and initial conditions. The grid is uniformly spaced with 51 points (i.e., 
NX = 51) over the length of the duct L. All results are calculated for an ambient base 
state in which P 0 = 101325, Co = 343, p Q = 1.2, and L = (throughout this work, 

the international units of measure axe used). The sound source oscillates sinusoidally at 
a fundamental frequency of ™ Hz (i.e., u> - 2200). Uniform time steps are used and 
the CFL number is ^ of its maximum value (i.e., cr t = 0.05). This time step is used 
because, in addition to satisfying the CFL criteria, it is small enough to provide accurate 
numerical data for evaluating the convolution integral in the radiation boundary condition 
(see equations (30) and (31)). The accuracy of the numerical solution at several points in the 
duct was investigated; and as expected, these results were more accurate than those results 
computed at the outflow boundary. Therefore, to limit the number of graphs, results are 
presented only at the outflow boundary x — L. 

The frequency- domain impedance Z(u> n ), which corresponds to outgoing plane waves, is 
the characteristic impedance of the medium for the current one-dimensional problem: 

Z(uj n ) = p Q co (38) 

Thus, in one space dimension, the temporal impedance (equation (32)) reduces to the product 
of the characteristic impedance of the medium poCo and a dirac delta function — r): 

Z(t - r) = p 0 co<7w(< - r) (39) 


a N (t - t) = - 



e »nu)(t-r) 


1 

T 


N 

1+2 ^2 cos tuv(t — r) 

n=l 


T 


sin - N * l u}(t — r) 


sin 




(40) 


Although the upper limit of summation N would be infinity, for a general wave the assump- 
tion of a band-limited periodic wave allows N to be truncated to a finite integer, which 
implies that the smallest acceptable value of N depends upon the bandwidth at the source. 
If the initial condition is transient free and the sound source is a pure tone, then a one term 
expansion (i.e., N = 1) of the Dirac delta function is sufficient. However, initial conditions 
with significant transients require a much larger value of N before the temporal impedance 
can be accurately represented. 

Note that because the temporal impedance Z(t — r) is a Dirac delta function in one 
space dimension, the integration in equation (30) could be performed analytically by using 
the shifting property of the Dirac delta function: 
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( 41 ) 


p(L,t) = p 0 cou(L,t ) 

Equation (41) can be rewritten in terms of p and u to yield 


P(L , t) - P Q = p 0 cou(L, t) 


(42) 


Equation (42) can be differentiated with respects to time to recover the linearized form of 
both the Hedstrom(reference [2]) and the one-dimensional Thompson (reference [3]) boundary 
condition: 


dp(L,t ) du(L,t ) 

~~dt '’ o0 ° Ft = 


(43) 


Unfortunately, in higher dimensional problems the temporal impedance will not be a Dirac 
delta function, and this integration must be carried out numerically. It is for this reason 
that sample calculations are presented using the exact analytical result (equation (41)) and 
results from the numerical integration equation (30). 


5.1 Results for pure-tone source 

In the first sample problem the source is a pure tone: 

g(t) = ,4sin(u>f) (44) 


A — ePo (45) 

where the parameter e must be considerably less than unity for the solution to exhibit linear 
behavior. Results presented in this section were calculated with e = 10~ 3 . The exact periodic 
acoustic solution, without reflections, can be computed from linearized acoustic theory: 


Pexact(x, t) = A sin (ut ) (46) 

Co 

This exact analytical solution is used as a test case to verify the accuracy of the radiation 
condition. We introduce the normalized error function E(x,t ) conveniently defined as the 
difference between the exact solution and the MacCormack solution: 


E>/_ 4\ Pexact(%it) PMacCormack{ x it ) / A rr\ 

E[x,t) = (47) 

The MacCormack solution contains error that is attributable to both discretization and 
computer roundoff as well as to initial transients that are reflected and not absorbed by 
the radiation condition. For a sufficiently fine grid, that is free of initial transients, the 
distribution of error E(x, t) will be considerably less than unity at the periodic state. Further, 
if the grid is sufficiently fine, then E(x, t) provides a measure of the initial transients that 
cannot pass through the radiation boundary. 

Graphs of the normalized outflow error E(L, t) for N = 2 and TV = 60 (see equation (40)) 
axe shown in figure 3. The curve that is labeled exact integration is obtained by using 
equation (41) as the radiation condition for MacCormack’s scheme. The initial condition 
that is used to generate each curve is transient free. The transient-free initial condition is 


9 



obtained by using t — 0 in equation (46). Because the initial data are transient free, E(x,t ) 
contains the effects of discretization and roundoff error only. The curve that corresponds 
to the smallest number of terms in the Dirac delta function expansion (i.e., N = 2) has 
the largest error. However, on this scale the effects of these errors axe small for the chosen 
grid since they axe less than 2 percent of the amplitude of the periodic state. Furthermore, 
when N = 60 the discretization and roundoff errors are just as small in magnitude as those 
obtained from the exact integration. Graphs of the normalized outflow acoustic pressure that 
is computed from the boundary conditions are compared with the exact periodic solution in 
figure 4. The solid curve denotes the exact periodic solution from linear theory; the broken 
lines denote the solution obtained with the temporal impedance radiation condition. The 
acoustic pressure has been normalized with the amplitude A and the physical time t with 
the fundamental period of the source T = Each curve collapses onto the curve for the 
exact periodic solution. 

The next example is presented to specifically show the effects of the radiation condition 
on initial transients. The result in figure 5 are similar to those in figure 3, when the system 
is initiated from rest (see equations (6)-(8)). This initial condition generates significant tran- 
sients, which must be eliminated from the solution domain by the radiation condition before 
a periodic state can be achieved. Because the discretization and roundoff error introduced by 
the current grid is small in comparison with the periodic amplitude (see figure 3), the error 
plotted in figure 5 can be attributed to transient effects. Approximately 0.8 of one period 
(i.e., t = 0.8T) is required before the initial transient reaches the outflow boundary. The 
effect of the number of terms N on the transient is apparent. When N = 2, the normalized 
error is unacceptable and grows in time. Thus, for a two-term expansion the transient wave 
does not pass rapidly out of the solution domain. However, when N = 60 the temporal 
impedance is more accurately represented, and a much larger portion of the initial transient 
is eliminated by the radiation condition than is eliminated when N = 2. This difference 
is especially true for large times. Figure 5 clearly shows that the best results are obtained 
when the time integration in the radiation condition is performed exactly. Furthermore, 
when exact integration is used the transient wave is transparent to the radiation boundary. 

Graphs of the outflow acoustic pressure which is computed by using the initial condi- 
tion with significant transients, are compared with the exact periodic solution in figure 6. 
Discrepancies between the exact and numerically computed curves in the figure are due to 
transient effects. Results calculated with N = 2 are not accurate. When N = 60, the 
initial transient is passed rapidly out of the solution domain by the temporal impedance, as 
expected. Superior results are obtained if the time integration in the radiation condition is 
performed exactly. After four time periods, the periodic solution obtained from each curve 
is accurate, except for the curve that corresponds to N = 2. 

5.2 Results for Broadband Source 

In the final sample problem the source is a multitone source: 

g(t) = A[sin(u>£) + sin(2tjf)] (48) 

where the amplitude A = 10 _3 Po and the fundamental circular frequency u is 2200 rad/s. 
Sources of this type are typical of those generated by ducted fans. The outflow acoustic 


10 



pressure is plotted in figure 7. The initial condition used to generate the results in this 
figure contains significant transients (i.e., the disturbance was initialized from rest), and a 
60-term expansion is used to represent the Dirac delta function. Results computed with 
the time-domain radiation condition agree extremely well with the exact results from linear 
theory. No evidence is found that indicates reflections of the initial transient on this scale. 


6 Conclusions 

A general theory for deriving time-domain nonlocal, nonreflecting boundary conditions for 
duct acoustics has been presented. The formulation is presented here in one space dimen- 
sion and time; however, this formulation can be extended in a straightforward manner to 
variable-area, higher dimensional, and acoustically treated ducts. The temporal operators 
for the time domain are shown to be convolution integrals of known frequency- domain oper- 
ators. Accuracy of the time-domain condition has been tested by applying a MacCormack 
integration to the Euler equations in one space dimension and time. The time limit for 
the time integration was arbitrarily set at 10 source periods to control computational cost. 
The MacCormack integration utilizes the boundary condition at an artificial boundary to 
limit the computational space and prevent spurious reflections. Results of this study show 
that the boundary conditions are extremely effective on both pure tones and broadband 
sources. Results obtained for a variety of initial conditions compared extremely well with 
known nonreflecting solutions after integrating the differential equations for no more than 
four source periods. Results show that a large number of terms are required in the temporal 
impedance operator when a significant amount of transients exist in the initial condition. 
When significant amounts of transients exist time steps that are smaller than the size estab- 
lished by the maximum CFL number must be taken in the numerical scheme. The smaller 
time steps are neccessary to accurately evaluate the convolution integral in the radiation 
boundary condition. 

The success of the temporal boundary conditions motivates their extension to higher 
dimensional and acoustically treated ducts. In higher dimensions, the boundary condition is 
more complex and involves a much greater computational effort. For example, the boundary 
condition will be a matrix convolution operator in two dimensions and a tensor convolution 
operator in three dimensions. 

Although a major purpose of the boundary condition is to improve accuracy and reduce 
the computational effort necessary for the solution of duct problems in large- volume domains, 
some point may be reached at which the added boundary-condition complexity, no longer 
results in a net reduction of computational effort when the computational volume is reduced. 
In this case, the major advantage of the boundary condition is its ability to account for the 
effects of sound-absorbing materials. 
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Figure 1: Infinitely long duct and coordinate system. 
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Figure 2: Continuous (x,t) and discrete domain ( Xj,t n ). 
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Numerical lntegration(N=2) 
Numerical lntegration(N=60) 
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Figure 3: Error function for pure-tone source (transient free). 
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Normalize outflow acoustic pressure, p(L,t)/A 


Exact steady state solution 
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Figure 4: Normalized acoustic pressure for pure-tone source (transient free). 

15 



Numerical Integration (N=2) 
Numerical Integration (N=60) 
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Figure 5: Error function for pure-tone source (significant transients). 
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Normalize outflow acoustic pressure, p(L,t)/A 
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Figure 6: Normalized acoustic pressure for pure-tone source (significant transients). 
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Normalized outflow acoustic pressure, p(L,t)/A 


Exact steady state solution 
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Figure 7: Normalized acoustic pressure for multitone source (significant transients). 
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